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Abstract 

We show that the reformulation of the de Rham-Gabadadze-Tolley mas- 
sive gravity theory using vielbeins leads to a very simple and covariant 
way to count constraints, and hence degrees of freedom. Our method 
singles out a subset of theories, in the de Rham-Gabadadze-Tolley fam- 
ily, where an extra constraint, needed to eliminate the Boulware Deser 
ghost, is easily seen to appear. As a side result, we also introduce a new 
method, different from the Stuckelberg trick, to extract kinetic terms 
for the polarizations propagating in addition to those of the massless 
graviton. 

1 Introduction 



During the past few years and in particular following discussions of the DGP model 
[T] and its cosmology [21 13] , there has been renewed interest in theories of "massive 
gravity" (see e.g. [HIS] for reviews). The unique consistent theory for a free massive 
spin-2 field was known for a long time to be the Fierz-Pauli theory [6]. This theory 
propagates 5 degrees of freedom of positive energy, one of which is a zero helicity po- 
larization responsible for the celebrated van-Dam Veltman Zakharov discontinuity: 
namely that, however small the graviton mass, Fierz-Pauli theory leads to different 
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physical predictions (such as hght bending) from those of hnearized General Relativ- 
ity Considering self-interactions of the massive graviton leads to a mechanism, 
first discussed by Vainshtein [S], which can actually restore the continuity towards 
well tested predictions of General Relativity |9l HOl [TTl [12] . However, massive gravi- 
ton self-interactions also introduce a new pathology: the fact that a ghost-like 6th 
degree of freedom propagates generically in the full non linear theory, as was first 
pointed out by Boulware and Deser ^J^. It was long thought impossible to obtain a 
massive gravity theory devoid of this ghost (see e.g. [121 HI])- However, a family of 
massive gravity theories was recently proposed by de Rham Gabadadze and Tolley 
(dRGT in the following) [151 UHl [13 in which the absence of ghost was first pointed 
out in the so-called decoupling limit [16] (using, in particular, the approach of Refs. 
[m [121 [IS]) and then fully confirmed at the nonlinear level by a Hamiltonian 
analysis, later extended to bimetric theories [221 [211 [221 [23] (see also [261 [27]). The 
Hamiltonian analysis of these models remains however complicated and does not 
clarify the reasons behind their soundness. 

Here, we will show that the reformulation of dRGT models in terms of vierbeins 
leads to a simple way of extracting covariant (Lagrangian) constraints. Some of the 
Hamiltonian properties of such a reformulation have already been analyzed in [2B] 
(building on the older work of [30j), where it was underlined in particular that one 
of the two necessary supplementary Hamiltonian constraints was easier to obtain 
with vierbeins than with the metric formulation. The analysis presented here is, 
however, different from the one presented in Ref. [28], in particular because it is 
fully Lagrangian and also because some of the arguments given there are completed. 

This paper is organized as follows. We first recall (section ^ how one can count 
degrees of freedom of massive gravity using covariant constraints. We then (section 
[3]) introduce the dRGT theories both in the metric and in the vierbein formulation. 
We then present our covariant way of obtaining constraints (section [1]), which is 
closely related to what can be done in quadratic Fierz-Pauli theory. As a side 
result, we also introduce a new trick, different from Stuckelberg's, to extract kinetic 
terms for the polarizations propagating in addition to those of the massless graviton. 
This trick, valid for a subset of vierbein formulated dRGT theories, is presented in 
appendix [A| 

2 The Boulware-Deser ghost from covariant con- 
straints counting. 

Let us first introduce the (quadratic) Fierz-Pauli theory [6J. This can be defined on 
a flat space-time by the following action (in the absence of matter source) for a 

^Note, however, that these results have been debated [24l 125) . 
6 See also [29]. 

^For simplicity, we only discuss in this section the case with D = 4 dimensions, but such a theory 
can be introduced in a similar way for arbitrary D, where it is found that a massive graviton in 
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rank-2 covariant tensor /i^j, 

Sh,m = -Ml J d^x [id^K,f - {d^hf + 2 {d^h) (d^K) - 2 {d^hj {d'^h^'p) 

[h^.h^'' - h^)] . (1) 

Here Mh is a mass parameter, all indices are moved up and down with a flat canonical 
Minkowski metric ?7^jy, and h is defined hj h = h^i^i]^^ . The terms on the right hand 
side of the first line of Eq. ([1]) are obtained by expanding the Einstein-Hilbert 
Lagrangian density y/^R{g) at quadratic order around a fiat metric, writing g^^, = 

+ ^Atf '^^^ mass terms appear in the second line on the right hand side of 
Eq. ([T]) and this particular combination of /i^ and {h^uY is the only one able to 
give a mass to the graviton h^y in a consistent and ghost-free way. Note that this 
theory explicitly breaks general covariance and also that it uses two rank-2 covariant 
tensors, h^i, as well as rj^^ which serves as a background on which h^^^ propagates. 

Fierz-Pauli theory can be non linearly completed by considering actions of the 
form 

Sg,^ = Ml j d^x^ [R{g) - mV (M)] , (2) 

where is a suitably chosen scalar function of Ai'^^ = g^" faui ^ and Mg are 
again mass parameters, R{g) is the Ricci scalar constructed from the metric g^j_v, 
and the theory contains, besides the dynamical metric g^^, a non dynamical metric 
ffj_,^ usually considered to be flat. If one wants to consider a (non linear) "massive 
gravity" the potential V should be chosen such that (i) when f^y is taken to be rj^^, 
d/Mu = ViJ.v is a solution of the field equations, and (ii) when expanded at quadratic 
order around this fiat background, the action ([2]) has the Fierz-Pauli form ([T]). Note 
that when one has two metrics, one can write any non trivial non derivative invariant 
built from the metrics as a function of A4, and hence the only restriction (besides 
diffeomorphism invariance) comes here from requirements (i) and (ii) on V. Note 
also that it is easy to figure out that there are infinitely many functions V that 
satisfy these requirements (see e.g. [31]). 

As first noticed by Boulware and Deser [13], quadratic Fierz-Pauli theory ([T]) 
and its non linear version ([2]) differ dramatically as far as the number of propagating 
degrees of freedom is concerned. Consider first the former theory. Varying action 
([1]) with respect to h^^, one easily obtains the field equations 

df,dyh + Uh^y - dpd^hP - dpd^hP + r]py{d''d''hp^ - dh) = {h^^ - hr]^^) . (3) 

Here the left hand side is just the linearized Einstein tensor and hence, as a conse- 
quence of Bianchi identities, its divergence vanishes. Thus, from the right hand 
side of Eq. ([3]) , we get 

V - duh = 0. (4) 

D dimensions, as described by the quadratic Fierz-Pauli theory, has {D — 2){D + l)/2 physical 
polarizations, i.e. the same number of polarizations as a massless graviton in D + 1 dimensions. 
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Taking another derivative of this equation yields 

d'^d^'h^^ -0/1 = 0, (5) 

where the combination appearing in the left hand side is just the linearization of the 
Ricci scalar. Thus contracting both sides of Eq. ([3]) with rj^'^ yields 

Uh - S^'S^ V = \'f^^h, (6) 

which together with (|5]) shows that h vanishes in vacuum. This in turn means, using 
dl]) that is transverse. The two equations we just obtained, namely, 

d'^K, = 0, (7) 
/i = 0, (8) 

together give 5 Lagrangian constraints (being first order) and this removes 5 of 
the a priori 10 dynamical degrees of freedom of /i^j,. Hence, quadratic Fierz-Pauli 
theory propagates 5 polarizations. A similar conclusion can also be reached using a 
Hamiltonian counting (see e.g. [T3]). 

In contrast, a generic non linear massive gravity propagates in addition a sixth 
ghost-like polarization. This was first argued by Boulware and Deser in Ref. [13] 
and thus the extra propagating mode is usually called a Boulware-Deser ghost. 
Schematically, this comes from the fact that the analog of the constraint ([8]) is lost, 
while there are still 4 constraints similar to ([7]). Indeed, now vary action ([2]) with 
respect to g^^ to obtain, 

G,, = mXl'f^ (9) 

where G^,y is the Einstein tensor built from the metric g, and T^^''^^ is obtained 
from varying the term coupling the two metrics \/—gV{A4), so that it contains no 
derivatives. Taking the covariant derivative (with respect to the metric g^^) of 
both sides of the above equation, gives 

V/^T^9'/) = (10) 

which therefore only contains first derivative^ (and in particular no derivatives of 
the metric f^u if this metric is just taken to have the canonical Minkowski form r]^;^) 
and hence yields 4 constraints on the dynamical metric g^j^^,. These constraints are 
the analog of ([7]) in the non linear case. On tracing over Eq.([n]) and using derivatives 
of Eq. flTU]) . there is now (as opposed to linear Fierz-Pauli theory) no reason to get an 
extra constraint (cf. also [E]). And in fact, it was thought impossible to construct 
a non linear Fierz-Pauli theory, with a suitable potential V, devoid of the Boulware 
Deser mode [T3l [T4] until the work of de Rham, Gabadadze and ToUey (henceforth 
dRGT) [151 [IE]. 



^Note that the derivatives of g only appear in the Christoffel symbols of g whenever a g covariant 
derivative hits an /^^ metric, since the action of such a derivative on g^^^ vanishes, and hence yields 
no derivatives. 
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3 The metric and frame formulations of de Rham- 
Gabadadze-Tolley theories 

3.1 Metric formulation 

dRGT theories are non-hnear Fierz-Pauh theories for which the function V takes a 
special form. We will use here the parametrization of dRGT theories proposed in 
Refs. [231 [32] • We begin by introducing the four functions Fi, F2, F3, F4 defined for 
an arbitrary n x n matri^j X'*^, representing elementary symmetric polynomials of 
the eigenvalues of X, and given by (see e.g. |23] ) 

FiiX) = [X] (11) 
F2{X) = \{[Xr-[X']) (12) 

FsiX) = ^{[Xf-3[X][X']+2[X']) (13) 

F,{X) = ^{[X]'-Q[Xnx'] + 3[Xr + S[X][X']-Q[X']) (14) 

where [X] denotes the trace X"^ of the matrix X. For general k, one defines Fk as 



Ffc(X) = ix«^[,,...X^^,], (15) 



where here and henceforth brackets [ ] denotes the sum over unnormalized antisym- 
metrized permutation j^. In particular, for a n x n matrix X one haj"1 that 

det(X) =F„(X). (16) 

Furthermore these functions appear in the expansion of the characteristic polynomial 
of the D X D matrix X^^,. Indeed, defining by 

n = x^. - (17) 

It follows that (here for D = 4) 
det (X'^^ - = F^iY) 

k=A 

= Y.{-xfF,_,{X) 

k=0 

= x^FoiX) - x'Fi{X) + x^F2iX) - xFsiX) + Fi{X) (18) 



^With a, a line index belonging to {1, n}, and 6, a column index belonging to {1, n}, and 
n having so far no relation to the space-time dimension D. 

""^^and similarly, parentheses ( ) will denote the sum over unnormalized permutations. 
^-^As a simple consequence of Cayley's theorem. 



5 



with the convention that Fq = 
by an action of the form ^23j 



1. The dRGT theory [151 IISI [H] can now be defined 



S 



Ml 



k=4 



k=0 



(19) 



where the /3„ are arbitrary parameters. Note that Pq just parametrizes a mere 
cosmological constant A, which by itself does not give any mass to the graviton, 
hence in the following we will generally trade /3o for A. It might nonetheless be 
necessary to keep a non vanishing Pq in order to have Minkowski space-time as a 
solution (and fullfill condition (i) of the previous section). On the other hand, the 
highest order term proportional to gives no contribution to the field equations 

of g^u, since (^Vd'^f) = V-9 det (^a/^-v) = a/^. Hence, in D = 

4 dimensions, there is a three parameter family of non trivial theories, indexed 
by parameters with n = 1,2,3. This can easily be extended to D arbitrary 
dimensions by considering actions of the form 



S 



dPx^—g 



k=D-l 



i?-2A-m2 l^kFkiyg^f 

k=l 



(20) 



where F^ are defined as in f lTBl) . Note that the above definitions (fTI?l) - fl2Ul) use a real 
matrix square root of the tensor Ai = g~^f- In general, however, there is no reason 
for this square root to exist for arbitrary metrics g^i, and /^,^ (see e.g. |37] and [S3]), 
and hence one has to assume that it does exist for the above definitions to make 
sense (we will come back to this issue later). When it does we define this square 
root as 7, and write 



(21) 



such that one has (for D = 4 dimensions) 



S = Ml (fx^^ 



k=3 



i?-2A-m2^/3fcFfc(7) 



k=l 



(22) 



An alternative formulation to action ( fT9l) is to consider 7 in f l22l) as an indepen- 
dant field and to enforce the relation ( 12T1) by a Lagrange multiplier c^*", adding to 
the Lagrangian a term of the form (in the spirit of e.g. Ref. [T5] ) 



(23) 



This alternative does not feature the presence of the unpleasant square root in the 
action. Whichever way is chosen (i.e. f ll9p or f l22p together with f l23p ) the presence 
of the square root directly in the action or via a Lagrange multiplier is a somewhat 
inelegant aspect of the theories considered. As we will show, the vierbein formulation 
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of these theories (or at least of a subset of them) offers a nice alternative which does 
not suffer from the same lack of elegance. 

For future reference, we also define actions by 



Mi 



Ml 



Mi 



(24) 



R - m^det [K^g-^f - 5^ 
j d^x^/^ [R - m^det (/t7^^ - ^^^J] 

-^[R-m" (Fo(7) - «:Fi(7) + '^'^^2(7) - '^'^^3(7) + '^'^4(7))] 



d x\ 



where k is a dimensionless parameter (and we have used identity ( |T8|) ). It is easy 
to see that on taking linear combinations of such models (with different values of 
k), one can obtain any model (fT9l) with arbitrary coefficients (in (using the non 
vanishing of a Vandermonde determinant). 



3.2 Vielbein formulation 
3.2.1 Generalities 

In order to formulate the dRGT theories in D dimensions using vielbeins, let us 
define E"^ and L"^ to be two bases of 1-forms obeying at every space-time point 

g^'E^^E^ = rL\L\ = , (25) 

or equivalently 

t^abE'^^E^, = g^, , (26) 

VabL^^.L'', = f^, . (27) 

We will also need the vectors ca and (.a, respectively dual to the 1-forms E^ and 
L"^, that verify 

E^(eB) = E^^CB^ = 5^B, (28) 
L'^iin) = L^'/b^ = (29) 

For future use, we define the (D—n)— forms E\^ ^ (using the notations of Ref. [35] ) 
by0 

EX...A^ ^ (D^'^-^-^^"" A ■ ■ ■ A , (30) 



-"^^Our convention is that Greek letters denote space-time indices, while capital Latin letters 
denote Lorentz indices that are moved up and down with the canonical Minkowski metric "qAB- 
We will, however, also use a particular coordinate system in which the /^i/ metric just takes 
the canonical Minkowski -q^i, form and the vierbeins L'^^ have components 'J'^^. In that case it 
sometimes turns out to be convenient to use the same type of letter to denote Lorentz and space- 
time indices - one then has to pay attention at the order of the indices to be able to discriminate 
between them. 

^■^Notice that the [D — n)— forms carry n Lorentz indices and that this definition also 

makes sense for n — in which case the form E* is just proportional to the volume Z?-form. 
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where eAi...Ao is the totally antisymmetric tensor verifying £123. ..d = 1- It then 
follows that 



n 



A = J2(-^r-'S\EX...A,^^A,,,...A. ■ (31) 

k=l 

Using the vielbeins E"^ the Einstein-Hilbert action for the dynamical metric g^^, 
reads 



M^'^ J Q^"" A EXb, (32) 



where 

^ ^^AB ^ ^A^ ^ ^CB ^ (33) 

is the curvature 2- form associated with the spin-connection . The latter is a one 
form taking values in the Lie Algebra of 5*0(1,3) (hence it is antisymmetric in its 
Lorentz indices A, B) and can be expressed in terms of the vielbein E"^ assuming 
(as we shall do henceforth) the torsion free condition 

VE"^ = dE^ + u^B AE^ = 0, (34) 

where the derivative operator V acting on an arbitrary p-form carrying Lorentz 
indices 11^^"'^'^^, o is defined by 

p=n 

p=i 

p=m 

- ^ C^^Bp A n^''"^"Bi...Bj,_iCBp+i-i?„ -(35) 



Using flMj) . one finds easily the components of the spin connection wabc which are 
defined by one of the two equivalent relations 



Wabc = e^^UABti, (36) 
3^ = E^ ^ 

and are given by 



UABf, = E'^^wabc, (37) 



1 

- eAVd^EBu - caVO^Ecu + eB^eA^d^Ec,), 



WABC =^{eB^ec''d^EAu - ec^CB^'d^EAu + ec^CA^d^EBu 



while the components of the curvature read 

e[c'^a^(w^^D]) + w^e\cw''''d\ - w^'^E'^.e^c'^d.eD]'' (39) 



AB \ I „,,A „,RB „,AB TpG 



i4\Ye recall that this condition can also be obtained by considering the spin connection as an 
independent field, and writing its field equations. 
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The curvature satisfies the Bianchi identities 

Pfi^^ = + uj^c A Qf'' + OJ^'c A r]^^ = . (40) 

The Einstein tensor has a simple expression in terms of the E\^(j forms and is given 
by the D — 1 form 

^ ^^^^^ A El^c- (41) 
Furthermore it obeys Bianchi identities 

VGa = , (42) 

which involve derivative of the D — 1 form Ga, and hence just yield D coordinate- 
scalar equations. Additionally, because of the local Lorentz invariance of the Einstein- 
Hilbert term, the Einstein tensor decomposed as Ga = Ga^E^ satisfies 

G[AB] = 0, (43) 

i.e. Gab is symmetric. 



3.2.2 Mass terms and field equations 

As discussed above the mass terms of the dRGT theory are expressed in terms of 
the (matrix) square root 7 of g~^f. Defining the Lorentz tensors ca^ (whose indices 
are moved up and down with tjab) as 



eA^ = eA' 



'L\, (44) 



a sufficient condition for this square root to exist is that the vierbeins obey the 
condition [251 131] (see also where it is shown that this condition is not necessary) 

e^^ = e^^, (45) 

in which case, 7 defined as 

I'u = e/L^. (46) 



verifies fl2T|) . Note that whenever the non dynamical metric f^y is flat, a convenient 
choice of vierbein can be made by first choosing cooordinates where /^j, takes 
the canonical form rj^y, i.e. 

ft^v = Vi^u, (47) 
and then choosing L"^ = dx"^, i.e. such that (in components) 

L\ = S\. (48) 
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In that case it is sometimes convenient to label space-time indices and Lorentz 
indices with the same set of letters, which we will do using latin capital letters. 
When the choice (H7ll - (H51) is made, one has 

eA^ = e^^ (49) 
and the constraint f H5]) simply reads 

= e^^, (50) 

stating that the vierbein e^'^ can be represented as a symmetric matrix. 

On substituting the expression (H^ into the mass terms of Eq. (HM it follows that 
these can be rewritten in terms of the vierbeins as 



n=0 



A---AL^" AEl . , (51) 



where we have absorbed irrelevant numerical coefficients n! by redefining the 
Also note that using the same substitution (and = det(£''^^) = £"), the action 
S'k of Eq. (12^ now reads 



= j (fxE [R - mMet [t^e/L^^ - 5^ J] 
j (fx [ER - mMet {kL^^ - E^J)] 



(52) 



= Mj, j (fx [ER - mMet (kL^^ - E\)\ . (53) 

As a side remark, we note that from action ( !53l) one can extract the kinetic terms 
for the extra polarization of a massive graviton in a very simple way i.e. by a simple 
shift of the vierbein. As far as we know, this trick has not been noticed before and 
differs from Stuckelberg's. However, since this is not the main subject of this paper, 
we discuss it in more detail in appendix \M 

Hence, generalizing to D dimensions, one is led to consider theories defined by 
the action 

D-l 

S = / A E*^B - M^-^m^ ^ /3„ / A • • ■ A L^" A , (54) 

where the kinetic Einstein-Hilbert term can easily be generalized to the Gauss- 
Bonnet-Lovelock terms in D dimensions. 

Varying this action with respect to the forms E"^ we get the following field 
equations in vacuum (see appendix [B] for a derivation) 

GA = tA, (55) 

with tA defined by 

1 

tA^-Yl Pn^^' A ■ ■ ■ A L^" A El^^,„^„ ^ tA^'E*^ , (56) 

n=0 
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and we have set to one for convenience. Using 

^ es^E^, (57) 
the coefficients can be computed to be 

D-l 

2 



n=0 



where 



3; ,B_Bi..._B„Ci...Ci3_„_i /'fin'\ 

~ (^D - n- 1)! eAAi...A„Ci...CD-„-i- loUj 

Note that the above set of theories (15^ (with field equations ( l55l) ) are perfectly 
well defined without imposing the constraint (1451) . Furthermore we will show that 
in some cases, the constraint fH5|) arises as a consequence of the field equations 
0. That this is the case has already been argued in Ref. [28]. However, the 
arguments presented in Ref. [28] use a decomposition of a generic vierbein (also 
adopted in Ref. [30]) as a product of a Lorentz transform with a symmetric vierbein. 
More specifically they assume the validity of the "Minkowski" version of the polar 
decomposition, assuming an arbitrary invertible matrix m can always be written as 
the product of a Lorentz transform A with a symmetric matrix s, 

m = As. (61) 

However, it can be shown (see e.g. [33]) that this decomposition only holds for a 
restricted set of matrices m, so that the arguments presented in [28] are in fact not 
fully conclusive. Notice that one has argued that the constraint 0451) can be set by a 
suitable Lorentz gauge choice [31] • In fact this argument also uses a decomposition 
such as fl6T|) . and so is not always valid. This will be discussed in detail elsewhere 
[33]. 

Hence, in the following we use action ( l54ll as a starting point, and begin by 
dropping the constraint ( H5l) . As such, the set of theories ( 15^ just defines some 
vielbein theory with a non dynamical vielbein and an unconstrainded dynamical 
vielbein E"^, carrying then, for D = 4, 16 a priori dynamical polarizations. 



4 Counting degrees of freedom using vielbeins 

In this section, we will successively see how various constraints reduce the number 
of physical polarizations of the a priori unconstrained vielbein. We stress that 

-"^^Notice that this dynamically enforces the existence of the square root 7, a nice feature of the 
vielbein formulation. 
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most of our results will be valid for any D space-time dimensions. For = 4, the 
constraints discussed in the following two subsections allow respectively to go from 
16 to 10 dynamical components and then from 10 to 6. We will show in the third 
subsection below that an extra scalar constraint can be obtained for a subset of 
theories. 

4.1 Constraints arising from local Lorentz invariance 

The constraints arising from local Lorentz invariance are encoded in the symmetry 
of Gab (see Eq. (143!) ) and read, from the field equations (155|) 

tiAB] = . (62) 

There are D{D — l)/2 independent constraints arising from ( l62|) . One therefore 
expects those to restrict the number of independent components of the dynamical 
vielbein e^'^, a priori D x D (i.e. 16 for D = 4 dimensions), to — D{D — l)/2 = 
D{D + 1) /2 which is the same number of components as in the metric g^^ (i.e. 10 for 
-D = 4 dimensions). We will even show that in some cases, constraints (!62|) further 
simplify, turning out to be equivalent to the conditions 

First, it is easy to see that assuming e^^ symmetric is sufficient to yield a 
symmetric t^^ (or equivalently Iab)- However, the converse is not true in general. 
For example the term proportional to ^2 in the right hand side of Eq. fl58|) yields a 
symmetric contribution to t"^^ when e"^^ is assumed antisymmetric. Of course, Eq. 
( H5l) can always be imposed (if needed) by introducing it by hand in the theory, i.e. 
by suitable Lagrange multipliers. As we just stressed, this will be compatible with 
field equations. 

There are, however, cases in which the equivalence can be established. In par- 
ticular when only /3d-i (besides possibly /3o, since the contribution of the latter 
in ( 15 8 p yields always a symmetric t^^) one has 

tA^ (X e^/^ . . . eB,_/-^C::S:.^ = det{L)E% , (63) 

where the matrix E is just the transposed inverse of e and is defined by its coefficients 
A satisfying 

£a = E^AeB. (64) 
Notice that with the choices ( 147|) - (H8|1 one has 

E^A = (■A^E'', = S^aE"", = E^'a . (65) 

In this case (non vanishing Pd-i), the symmetry of t^B implies the symmetry of 
Eab (because det(L) 7^ 0). This automatically implies that e^^ is also symmetric. 
This also happens when only /3i 7^ 0, in this case one has 

tA^" oc es/^^fS = tr(e)5f - , (66) 
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which again yields a symmetric e^^ . For an arbitrary combination of mass terms 
with non vanishing /9„ no definite conclusion can be drawn, however choosing suitable 
Pn results in the mass term appearing in Eq. and yields (with the choices (H7|l - 

)) 



tA^ oc ^"^"gg^^ ^"^^ = det{Ke - Id)[(/.e - Id)-^^ , (67) 



which includes the two previous cases as limiting cases (when k = and k — )■ oo). 

Notice that whenever we choose the gauge (HHI) . we have e^^ = e^^ and Eab = 
Eab, and thus the symmetry of the hatted quantities implies the symmetry of the 
vielbeins themselves. 

4.2 Constraints arising from difFeomorphism invariance 

Using the Bianchi identity ( H2|) in the field equations we get that 

VtA = 0, (68) 

which reads explicitly 

D-l 

2 



1 

= o 5Z ^^rrVL^' A • ■ ■ A L^" A E\a,...a^ = 0. (69) 



n=l 



Here, and for the remaining of this subsection and the next one, we choose the gauge 
^8|) which results in the vanishing of dL"^, the above equation becomes 



D-l 



J2 nPn io^'B A A L^^ A ■ ■ ■ A L^" A E*aa,...a„ = 0. (70) 

n=l 

On using L^' = es/' E^' = es/' E^', u^^b = w^^bc E'^ and the identity 
(deduced from (160|1 ) 

E^^A-.-AE^-A E^^...^,^ = e6^l:±, (71) 

where e denotes the volume D-form E^ A E^ A ■ ■ ■ A E^, we can rewrite fITU]) in the 
useful form 

D-l 

J2 w^^B cCB^'^eB/' • • • eB/"5S.:t ^ = 0' (72) 

n=l 

Constraints (168!) clearly remove another D degrees of freedom. So, in D = 4 di- 
mensions we have so far shown that from the 16 components only 16 — 6 — 4 = 6 
are left dynamical and one needs one more constraint to have only the sought for 5 
dynamical degrees of freedom of a massive graviton. In the following, specializing to 
D = 4 dimensions, we will discuss how such a constraint can arise in some specific 
cases. 
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4.3 Extra scalar constraint 



In this section we particularize our discussion to D = A dimensions and to the gauge 
choice fHHj) . The way we proceed to obtain the extra constraint is very similar to 
what was done in the quadratic Fiez-Pauli theory, as explained in section [21 where 
the extra constraint is given by Eq. (|8]). Namely one uses the constraints coming 
from the Bianchi identity (ISSj) (cf. Eq. (jl])) back into a suitable trace of the field 
equations fl55]) (cf. Eq. Q), which in our case reads 

A = A t^, (73) 

where m"^ is a suitable collection of one forms (labelled by the Lorentz index A). 
As we will now show, however, things proceed differently depending on the values 
of the (3n coefficients. 



4.3.1 Only /3i ^ 

We first discuss the case in which the only non vanishing coefficient /3„ is /3i (and 
possibly /So parametrizing then a non vanishing cosmological constant), because this 
case is the simplest. Then, we can rewrite the Bianchi identity fl72l) as (the term 
proportional to /3o automatically vanishes) 



(74) 



which immediately implies 



w'^AB = 0, (75) 

yielding then 

eB^'dcE'^A - CB^'dAE^'c = 0. (76) 
Choosing = E^ and tracing over the field equations fl55l) we get 

E^ AGa = E^ A tA. (77) 
The left hand side of the above equation can be rewritten using (jSHD and (|7T]) as 

E^AGa = ^E^An""^ AE\^c, 

~ £e^^9^(w7^^B) , (78) 

where the symbol ~ means that we only write terms containing second order deriva- 
tives and omit an overall constant factor. But the constraint (176]) tells us that these 
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terms vanish, so equation ( !77|) contains only derivatives of order one at most, and 
hence represents an additional constraint. This constraint, together with the four 
vector constraints (175]) can be elegantly rewritten using the following decomposition 
of the vierbein (valid only on the "branch" of vierbeins with positive determinants) 

= e^'E^ (79) 

with det(-E') = 1, which immediately implies 

E^B = e'^E'^B , (80) 
= e-'^CB^. (81) 

Now the vector constraint (176|) becomes 

Oact = ^eB^9[c^^A] (82) 

which, using the fact that the matrix whose matrix elements are the E^b has a unit 
determinant (which implies that cb'^OaE^c vanishes) can be rewritten as 

dA(T = ^CB^'dcE^A ■ (83) 
Eq. (1791) implies on the other hand that 

e^WABC = wabc - e[A^df,a t]b]c , (84) 

where the coefficients wabc are defined as those, wabc, of uab (see Eq.f l37|) ) by 
Lj'^B = w^BcE^ ■ In terms of those coefficients, the vector constraints ( 176|) just read 

SeA^^^BfT = -dBCA'' = -W^AB. (85) 

On the other hand, using (1791) and (18^ the scalar constraint (1771) just reads 

- WABcw^""^ + Iw^'abw^'^c + 2w^^ceA''eB''dHE''G = 2/3oe''^ + ^/3ienr(e^^) 

(86) 

which shows clearly that beyond the first derivative already constrained by the vector 
constraints another independent Lorentz invariant relation between first derivatives 
can be obtained. 

4.3.2 Only (32^0 

The starting point here is the same as above, starting from the vaninshing of T>tA, 
and using first only the symmetry of t^^ which implies 

w^BctA'' = 0, (87) 
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we can rewrite Eq. fl72|) . as 



= VtA = (32W^'BceB,''eB,^'5'^^l^le (88) 
-w'^cBeA^'en'' - w"" cDeB"" tA^) e (89) 

yielding 

w^CBCij^eA^ + w'^cDeA^CB'' - w'^cBeA'^eD'' - w'^cDeB'^eA'' = , (90) 
which reduces to the simpler form 

w'^cBeA^ + w'^AceB'^ - w'^ABec^ - w^'cAeB'^ = 0. (91) 

Due to the form of this equation and to the fact that the left hand side of flSSl) does 
not contain any it appears natural to use here = to trace over the field 
equations. Using the same logic (and notation) as before we get that 

L^AGa ~ eeA^d^iw^BC^'^' + w^^ceB'' -w^^Bec"") 

~ £eA^a^(w^c^eB^). (92) 

The right hand side above does not vanish in general, and hence, in this case the 
equation A Ga = A Ia does not automatically provide the extra constraint 
we need. However, if one imposes e^^ to be symmetric (which, as we stressed in 
14. H does not contradict field equations), then w^c^cb'" vanishes (by the same logic 
as in Eq. ( 187|) ). and hence we obtain an extra constraint similar to the one found 
above. However it should be stressed that this constraint has been obtained using 
a procedure which differs in its details from the one of 14.3.11 first we had to trace 
with a different set of one forms and second, we had to impose by hand the 
symmetry of e^^. 



4.3.3 Other cases 

A scalar constraint for an arbitrary combination of terms with non vanishing /3o, Pi 
and P2 can easily be obtained along the previous lines. One just needs to trace over 
the field equations with an appropiate combination of and in order to make 
all second derivatives disappear. 

However things proceed quite differently whenever is non vanishing. Indeed, 
let us now discuss the case in which only fS^ differs from zero. First recall that 
now there is no need to assume that e^^ is symmetric since it appears as a mere 
consequence of the field equations. The vector constraints also take a simple form, 
indeed we can rewrite ( 17U]) as 

3 

= VtA = -PsOO^^Bt^dx^ ^dx^ Adx"^^ Adx^'eAAiA2A3 

oc Psw'^AcE^Be. (93) 
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Hence, the vanishing of T>tj^ yields that of acE^ b which we can rewrite in terms 
of the vierbeins as 

dA [E^'b) - ObE^a = 0. (94) 
However in this case, it is not possible (as shown in appendix [U]) to find a collection 
of one forms m"^ such that the equation m"^ A Ga = A provides an additional 
constraint in the same way as in the previous cases (at least under some fairly general 
hypotheses on m^). Hence this case does not appear as transparent as the others 
as far as the existence of the extra constraint is concerned. The same would be true 
for an arbitrary linear combination of mass terms where ^3 is non vanishing. 



4.4 Recovering quadratic Fierz-Pauli theory 

It is instructive to show how the constraints of the quadratic Fierz-Pauli theory, 
Eqs. dZD-dH]) can be recovered from the constraints derived in sections 14.21 and 14.31 
by expanding the non linear massive gravity around a flat space-time. This will in 
turn allow us to show that the constraints we obtained are independant from each- 
other, as they should be. To do so, we first need to make sure that flat space-time is 
indeed a solution of massive gravity and provides a suitable background which might 
require adding a non vanishing cosmological constant. This addition does however 
not change any of the conclusion of the previous sections, as we already stressed. 

Let us then look at the linearized limit of the constraints. We expand the dy- 
namical vierbein as 

= dx^ + E{[), (95) 

CA = dA + e^A^ (96) 

where here and henceforth an index (i) denotes a first order perturbation. Writing 
the metric perturbation as h^^,, which verifies g^^ = r]^^ + h^^, we have that 

-^(^B) = ^AB, (97) 

e|f)^^ =-h^^ (98) 
Et,)B =-eg)^ (99) 

where, because we are working at linear order, Lorentz and space-time indices are 
identified. Choosing then as in (HHl) . and writing 

CB^ = 5^B + 4'^^, (100) 

we get that 

= ^i)"^ = -^^"^ • (101) 

Now at linear order reads t^^^ oc Tr(e(i))(5^yi — e^"*^ ' which implies that e^^ 
is symmetric and therefore that 

E^l = ^ . (102) 
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In order to be able to write the constraints in terms of h^B we also need the expres- 
sion of the connection 1-form. From f l38p we get 

''J^ABC = ^{dBhAC - dAhsc) ■ (103) 
It is now easy to see that, at linear order, the constraints ( !72|) are just 

wf,-^AB = \{dAh''B-d''hAB)=Q , (104) 

which is nothing but the vector Pauli-Fierz constraint (jlj. 

Let us now examine the additional scalar constraint corresponding to the scalar 
Fierz-Pauli constraint (j8]). To do this we look at the trace (|73ll . where m"^ is either 

or L"^ (or a combination thereof). As we would expect, at the linear level, the 
left hand side of Eq. ( |73l) vanishes 

A oc dA{wf^fB)e = ^d^idAh"" b - d''hAB)e = 0; (105) 

as for the right hand side, if the 1-forms m"^ are just L"^ or E"^ (or any combination 
of these) then one has at linear order 

A oc tr(e(i))e = h^A^^, (106) 

which yields h'^A = 0, as expected. 

5 Conclusions 

In this paper we have analyzed a vielbein formulation of a family of ghost free mas- 
sive gravity theories obtained from the de Rham-Gabadaze-Tolley metric theories. 
Summarizing here only our results for D = A dimensions and a non dynamical metric 
f^j^jj which is flat, this family is a three parameter set of theories characterized by an 
arbitrary combination of three mass terms (l^Tj) each proportional to a constant /3„, 
n = 1, 2, 3, in addition to a possibly non vanishing cosmological constant A (which is 
also given by a similar term proportional to /3o). In the vierbein formulation, start- 
ing from an arbitrary vierbein (hence with a priori 16 free polarizations) we showed 
that the constraints associated with local Lorentz invariance (encoded in the sym- 
metry of the energy momentum tensor) and those (vector constraints) associated 
with the Bianchi identities generally reduce the number of physical polarizations 
to 6. We then studied the three cases in which all but one among (32 and /^s 
vanish. In the first case (only /3i 7^ 0), we began by showing that the field equations 
impose the vierbein to be symmetric; and then, we showed how an extra scalar con- 
straint (a scalar combination of the field equations only containing first derivative) 
can be obtained, analogous to the tracelessness of the graviton field in the standard 
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(quadratic) Fierz-Pauli theory. This is obtained by taking into account the vector 
constraints as well as a suitable trace over the field equations. This extra constraint 
further reduces the number of propagating degrees of freedom in agreement with 
results obtained using the full Hamitonian analysis of the theory |20l |2T1 [221 123] . 
In the second case (only P2 non vanishing) a similar constraint can be obtained, 
with however the following two important differences: firstly one has to impose the 
symmetry of the vierbein which is not garanteed anymore by the field equations, 
secondly one uses a different trace. In the last case (only non vanishing), we were 
able to show that, even though the field equations do impose the symmetry of the 
vierbein, the previously followed procedure to obtain a scalar constraint does not 
work. The same would hence be true for an arbitrary combination of the three mass 
terms as soon as does not vanish, however a scalar constraint can be obtained for 
such a combination and a vanishing P3. Let us stress that is it quite remarkable and 
non trivial that such a scalar constraint can be obtained at all in some fairly general 
cases. However, the puzzling difference between cases does not by itself invalidate 
the results of [20] [2T| [22| [23] which are claimed to be valid for an arbitrary theory 
in the dRGT family (even though it might also open a way to reconciliate the con- 
tradictory claims which have been made in the literature about the Boulware Deser 
ghost). Indeed, it might just be that the extra constraint cannot be written in the 
most general simple space-time scalar equation involving only first deriva- 

tives (which is a stronger requirement than the existence of an extra constraint on a 
subset of propagating degrees of freedom). This requires more work and could also 
be checked by a proper Hamitonian analysis of the vierbein formulation. 
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A A simple way of extracting propagating de- 
grees of freedom with vierbeins 

Our starting point here is the subset of theories with actions (in vacuum) given by 
5*^ of Eq. fl53|) . It is well known that the Einstein Hilbert piece of this action can 
simply be rewritten as a term quadratic in the first derivatives of the vierbeins, such 
that the 5*^ action is also given by 

= Mj, j d^x[E G^''^'''^PEA^,uEBx,p-m'det{KL^,-E\)] (107) 

where G^^a'^^V 

is a (mixed: space-time and Lorentz) tensor containing vierbeins 
E^ (and their inverse) but not their derivatives. Hence, when one considers a fixed 
background vierbein L^, one can redefine the dynamical vierbein to be 

E^^ = E^^ - kL^^, (108) 

so that the above action becomes 

= Mj,j rfS [e G^^^^"''Ea^,,Eba,p - mMet (e^,)] + 5™, (109) 

where Smix involves a mixing between derivatives of E^ and (as well as a term 
containing only derivatives of L^), and G^B^vXp obtained from Q^B^vXp -\^^ replac- 
ing E^ with E^. When = dx^, Smix is simply given by 

Smix = Ml j d^xEK^^^^'^PEAp^uEBKp (110) 

where i^^^/^'^V depends on E^ and but not on their derivatives. Notice now 
that the first two terms in f ll09p are invariant under diffeomorphisms (the mass 
term has been replaced by a cosmological constant), and that the term f lllOp can 
just be interpreted as a modification to the kinetic term of the massless graviton, in 
particular it explicitly depends on the background vierbein L^. As such (together 
with the first term in f llOQp it encodes for the kinetic terms of the extra propagating 
polarization of the massive graviton. 



B Field equations in the vielbein formulation 

We vary the action fl5^ with respect to the vielbein E^. Note that for simplicity 
we have set Mp = m? = 1. This yields 

SS= f 5Vl^^ AEIb+ j A5E\B-^Pn f L^^A---AL^"A5E:^^...^,^ . (Ill) 
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Using that 5E\ = A ^aAi...a„ the last term can be rewritten as 



;il2) 



n=0 

The same method shows that the second term reads 



j SE^ A A E\j^c = 2 j SE^AGa. (113) 

Thus the only potential complications arise from the first term. Let us look at it 
more closely. Using the definition of the curvature 2-form as well as an integration 
by parts we can rewrite it 

J 5uj^^ A dE\B + j ^u^c A UJ^^ AE\b + j Lo^c A Soj^^ A E^^ ■ (114) 
Using the following straightforward consequence of the torsion-free condition 

n 

= -T.^^'a, a i?l„„^,_,BA,...A„ = 0, (115) 

k=l 

we see that 

dE\^ = dE"" A E^BC = -^^D A A E\^c = a A E*^^ + ou^ b A E^^ (116) 
and the expression (11141) becomes 



du"^^ A u^A A E*cB + / A tJ^B A E^^c 



+ / Su^c A o;^^ A + / A (5a;^^ A E* 



(117) 

which is clearly equal to zero. Therefore we have proven that 

6S = 2 J 5E^ A{GA-tA) , (118) 
which naturally implies the field equations f l55|) . 

C On the case ^ 

We show here that when is the only non-vanishing /3 (cf. subsection I4.3.3P 
there is no family of one forms which are polynomial in E^ and such that 
tji^AGa = m^AtA provides an additional constraint in the sense discussed above, i.e. 
such that A Ga does not contain second order derivatives. Using = m^sE^ 
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and the ~ notation introduced earlier we have, for an arbitrary collection of one 
forms 

m^AGA = A Q^^ A E*^^^ 

= ^m^nW'^'EFE'' AE'^AE^A E^^^ 

= ^m\W^^EF6^E5e (119) 

~ en''d,im^Ew'''^F)6'^EE^ 

~ eA^d^im'^Bw''^c + m^'cw^^B - nf cw^"" B)e . 

We show then that it is not possible to find (provided it is only built from 
and E^) such that the term in factor of /Sse in Eq. (p3|l can be written as the term in 
the parentheses above. More specifically, we prove that there do not exist invertible 
matrices Ai = {m'^B)i<A,B<d and S = {s'^B)i<A,B<d such that for all E^b 

s cw dJ^ b — ^ bw c + ^ cw b — b- \^^^) 

Indeed, let us first rewrite the above equation in the more convenient way 

{s^cE^'b - m^B^^'c + m^B^^c - m"^ e5^ c^"" b)w^'' d = 0, (121) 

and notice that we can treat w^^ d as an indeterminate (this because one can vary 
the derivatives of the E^ without varying the E^ or the e^). The first subtlety that 
we have to take into account is that w^^ o is antisymmetric in its first two indices. 
A second subtlety arises as a consequence of the symmetry of e^^ or equivalently 
of its inverse. Indeed in this case w^^ d acquires a new symmetry which reads 
W[B\cD\E^ eE^ F] = (this relation is just a consequence of the differentiation of 
TjAsE^ A dx^ = 0). Therefore the identification reads 

s^[ci?^B] -m^[^5^c] +m^[^5^c] -m^^^^^^^B] = VG[BVc]HA^Ke''''^'^E'' eE^'e , 

(122) 

where A^^^ is an arbitrary rank two tensor which is a function of the E'^b- Taking 
two traces over the indices A, B and C, D simultaneously and introducing the 
matrix £ = {E^B)i<A,B<d we get 

tT(M) = -(tr(5)tr(£:) - tr(^5)) . (123) 
6 

Now, taking only one trace first over the indices A, B and then over C, D, and 
plugging in the above result we get 

M = -{£S - £ti{S)) + -(tr(5)tr(£) - ti(S£))I - -X 

2 6 2 (^^24) 

= l-iS£- StTi£)) + ^(tr(5)tr(f ) - tr(5^))/ , 
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where X is the matrix whose matrix elements are the 

flG[AflC]H^^ K^^^^^ eE'^ F = ^GK^^^^^ EceE^ F — VGC^HK^^^^^ E^ eE^ p ■ 

(125) 

But since the A"^^ are functions of the E^b then necessarily the terms ^gk and 
E-^HK niust be symmetric in their free indices. Thus the right hand side of Eq. 
( I125p must be identically zero which shows that we may disregard the additional X 
term in fll24p and obtain 

SS - StT{S) = SS- Sti{S) . (126) 

The only S which verifies this relation whatever the S is S = aS. This gives us an 
associated matrix Ai and it suffices to verify that this combination does not work 
in order to complete the proof. 
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